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This is an English translation of my 1977 Russian preprint. It contains the first explicit definition 
of the pion distribution amplitude (DA), the expression for the pion form factor asymptotics in 
terms of the pion DA, and formulates the pQCD parton picture for hard exclusive processes. 
Abstract of the original paper: 

The large Q 2 behavior of the pion electromagnetic form factor is explicitly calculated in the non- 
Abelian gauge theory to demonstrate a field-theoretical approach to the deep elastic processes of 
composite particles. The approach is equivalent to a new type of parton model. 



At present {1}, it is usually accepted that the hadrons are complex objects composed of quarks interacting 
through a mediating gluonic field {2}. The approximate scaling in deep inelastic scattering indicates that the 
effective coupling of the quark-gluon interaction is small in the region of large transferred momenta. Hence, there 
is some justification for the hope that the quantum-field calculations based on the use of perturbation theory can 
be appied in the region of large momentum transfers. An efficient technique based on the use of the operator 
expansions and the renormalization group (RG) methods was applied to the study of the deep inelastic scattering 
processes [1]. Incorporating the recent studies of the pion form factor behavior at large Q 2 with the help of the 
operator expansions [2] and the analysis of the Feynman graphs in the a-representation [3], we investigate deep 
elastic processes involving composite particles with the help of a technique analogous to that used in [1] . 

Consider the pion electromagnetic form factor, treating the pion as a bound state of spin 1/2 quarks. The 
form factor of such a system can be written in terms of the Bcthe-Salpeter wave function \ as shown in Fig. 
la. Due to property \ = K ® X; where \ is the Bethe-Salpeter kernel, we can rewrite F v as shown in Fig. lb. 




b) 

Fig. 1. Electromagnetic form factor of a composite pion. 



Each graph can be formally written in the a-representation (see, e.g., [4], [5]): 

R-(P,n- S '(^)"""f D^ G »(«,P,P')exp{ i ,^ + i /( C »m;)} , (1) 

where the functions D,G,A and / are determined by the topological structure of the diagram. Note that X:X 
vertices should be treated as subgraphs rather than points. One of the most important results of the analysis of 
the Feynman graphs in the a-representation is that the large-Q 2 behavior in a theory with dimensionless coupling 
constant is determined by integration over a region of the a-parameters belonging to subgraphs the contraction of 
which into point eliminates the dependence of the graph on the large variable Q 2 . For the pion electromagnetic form 
factor, integration over small Ay (Ay = a a ) gives the leading contribution Fy ad (Q 2 ) ~ Q~ sti+2 , where ti is 

<t£V 

the twist (= dimension minus spin) of the field describing the i-th external line of the subgraph V, and Q = \f~-q 2 ■ 
Taking into account that, in a theory with vector gluons, t = 1 for quark fields, but t = for vector fields (both for 
gluons and the photon), one may conclude that the leading contribution can be given only by subgraphs having 4 
external quark lines and an arbitrary number of gluon lines (Fig. 2). 
Using the Mellin representation 

FAQ 2 ) = 77- / HJ){Q 2 ) J d.j , (2) 

one can say that integration over the region < Ay < 1/fi 2 gives a pole (J + l) _1 (l//i 2 ) J+1 / ®E®f, where E is 
the contribution due to the subgraph V and /, / are the contributions of the left and right weakly connected parts 
resulting from the contraction of V into point. 
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Fig. 2. Subgraph producing the leading large-Q 2 contribution. 



The functions /, / do not depend on J because the contraction of V into point eliminates the dependence 
on Q 2 . Summation over all possible subgraphs looks like a perturbative series, with E(Q 2 / fi 2 , g(fi 2 )) describing 
short-distance interactions while / and / correspond to the properly normalized quark-pion vertices: 



F ^Q 2 ) = ^E E ffa\g,P')E!? n (QVf,9)tf{v<',9,P) 



(3) 



n— m— 



where /, K denote the spinor indices. The leading contribution can be obtained only when / and / are projected 
on the axial structure [2,3,6]. The functions f£ = a n are defined through matrix elements of the twist- 2 operators 



2"*"- 1 (0|^(0) 75 {7m 1 £>m 2 • ■ • ^ n+1 }V(0)|P) = • • • P„ n+1 } a n (» 2 ,g(» 2 )) + \ ' 



(4) 



where is the covariant derivative acting on the quark fields, { } means that the corresponding function 
is symmetric in fii,.. . , /i n +i an d all its contractions with g^to are zero, and ^ 2 serves as the renormalization 
parameter. 

Applying \idjd\x to both sides of Eq.(4), we obtain a renormalization group equation 



oo oo 



n—n' m—m' 



m'm^nn' 



(5) 



where 



^ + ^dg 



a n (p 2 ,g(fi 2 )) = ^ z nn>a n > ■ 



(6) 



According to (6), there is mixing between operators with different spin (but the same twist). 

Let us illustrate the method of calculations on the example of the lowest-order diagram (Fig. 3a). The a- 
representation for the leading contribution of such a diagram in a non-Abelian gauge theory has the following 
structure: 



/ 2 \ Nl+n r r 



Vl,Vr 



D 2 (V L ) 4 W7 " l ' i) j £)2(V^ 



(it) 



; ^ (fi) >"4)isp{ 7 57AG fl } 



I dad/3 Sp {75 7 "7 Q 757 A 7^7n e ^^+/3 W+°(A 2 ) 



(7) 



where a4^ € V^a^ <E Vr. We also have taken into account that, for the leading contribution, the function Gl 
(G r ) depends only on P (P 1 ) and that D = D(V L )D(V R ) + 0(A), where A = a + [5. The meaning of the 2-trees 
L,L',R,R' is clear from Fig. 3b. 
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a) b) 

Fig. 3. Structure of the lowest-order diagram. 



Integration in the ncighboorhod of a <~ 0, j3 <~ gives 



v L ,v R 



D 2 (V L ) \-L-/D{V l ) 



G A (V L ) 



D 2 (V R ) 1-R-/D(V r ) 



G A (V R ). (8) 



We used the fact that for the diagrams of Fig. 3b type R + = D(Vr), where R± = R' ± R, and hence R/D(V R ) = 
(1 - R-/D(V R ))/2. Expanding into series 1/(1 - #_/£>(Vr)) = £(#_/£>)" wc obtain an expression of Eq. (3) 
type. The factor (R-/D) n just means that the corresponding operator contains n derivatives. Due to symmetry 
between R and R' , the contributions corresponding to odd n vanish. Let us introduce the partonic wave function 



Since \R_/D\ < 1, the function vanishes for |£| > 1. 

In other words, the function £>(£,^ 2 ) is such a function, the n-th moment of which is equal to a„(/i 2 ) {3}: 



(9) 



1 i _i_ ("—El™ 



(10) 



Hence, <,?(£) = ¥>(-£)■ 

The function 99 can be normalized by noticing that the magnitude of the matrix element of the axial current is 
known from the tt — > \iv process: 



(11) 



with /jr w 100 MeV. Since the axial current has zero anomalous dimension, the relation (11) is valid for any /i 2 . 
That is why we take £5(£,/z 2 ) = ^f^a^,^ 2 ). The final result then has the following form 



K (Q ) = ^™ S (Q 2 



ft C 2 {R) 
Q 2 N c 



2\|2 



(12) 



where 



1 . 



The representation (3) is equivalent to the following parton picture for deep elastic scattering. The splitting of 
the pion into two quarks with momenta x\P, x 2 P is described by the wave function <p(xi, x 2 , [i 2 ) 5(1 — x\ — x 2 ), while 
the fusion of the quarks with momenta y\P' , y 2 P' into the pion is described by the wave function tp*(yi,y 2 , fi 2 ) 5(1 — 
yi — 2/2)- The amplitude for xiP,x 2 P — ► y\P' ,y 2 P' transition is constructed according to the usual rules: 



FAQ 



">-/ 

Jo 



dx dy (p*(y, 1-y, y 2 ) E(x, y, P, P' , fi 2 ) <p(x, 1 - x, (i 2 ) 



(13) 
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For the diagram 3a {4} 



xyq z xq z 14 4 j 

the wave functions are symmetric ip(xi,x 2 ) = ip(x 2 ,x\) and — ip (j-^i ~^^) ' ^his P n y s i ca l interpretation 

suggests that ip(^,y?) is maximal for £ = and vanishes for £ = 1. Taking for simplicity |a(£,/i 2 )| = |(1 — £ 2 ) 
{5} (where 3/2 is determined by the normalization condition) and incorporating that the analysis of scaling 
violation in deep inelastic scattering under the assumption that the asymptotic freedom takes place, i.e., that 
a s (Q 2 ) = 9i n (Q2/ A 2) , gives A = 0.5 GeV [7-9], we obtain that a s (Q 2 = 2GcV 2 ) = 0.7, and, taking into account 
that C 2 (R) = 4/3, N c = 3, we find i^(Q 2 = 2) = 0.18. The formula F n {Q 2 ) = (1 + Q 2 /0.47)- 1 (the extrapolation 
of the p-meson pole well desribing experimental data) gives F^(Q 2 = 2) = 0.19. The factor m 2 /Q 2 determining 
the magnitude of the pion mass corrections is small for Q 2 > 2. Namely, m 2 /Q 2 < 1% . Thus, the corrections of 

the ^-scaling type [10] may be neglected. There remains the contribution of the twist-3 operators V'Ts D ■ ■ ■ D ip- 
It can be calculated in the same way as the axial contribution. One should also include the next term in the 
expansion of the function E(l,g(Q 2 )) having the order of magnitude of a s (Q 2 )/TT ~ 20%. 

Thus, we observe that the agreement of quantum chromodynamics for a s (2) = 0.7 with experimental data is not 
bad, at least in the region 1 < Q 2 < 3 GeV 2 . 

The parton formulation of our approach can be easily applied to other hadrons. The proton in deep elas- 
tic scattering is described by the wave function ip(x\, x 2 , x 3 , fJ 2 )S(l — x x — x 2 — x 3 ) related to operators of the 

ipiC'~f5'~ft_ l ip2{D23) N2 (Di3) Nl ip3 type, where C is the charge conjugation matrix. Unfortunately, there is no a priori 
information about normalization of the function ip(x±, x 2 , x 3 , /i 2 ). It can be extracted from the data on the proton 
form factor, and then the parton wave functions can be used to fit the data on deep elastic wide-angle hadron 
scattering: 

M(s,t)\ AB ^ C n = J \{dx^dx^dxWdx^5 (l - $>£ A) ) « " E^) 5 { l ~ E*£°) 5 { l ~ E^) 
x v^x^) V ^x^)e(x^\x\ B \x£\xI d \P a ,P b ,P c ,P d ^ 2 ) V {x[ A) )M B) ) ■ (14) 

The lowest-order approximation for E reproduces the quark counting rules [11,12] da/dt \ab^cd ~ s~ N+2 , where 
N = ha + Us + nc + no- In addition, there appears a factor [a 8 (ut/s)] N ~ 2 that gives an extra dependence on 
s. In the region ut/s ~ m 2 one should take into account the proton mass corrections. The function Q 2 F 2 (Q 2 ) 
behaves like a s (Q 2 ) even for moderately large Q 2 . For this reason, the investigation of the pion electromagnetic 
form factor seems to be the best tool for the experimental study of the nature of the quark-gluon coupling constant 
renormalization. The wave functions ip(x,Q 2 ) also depend on Q 2 due to the anomalous dimensions. However, 
preliminary estimates show that the change of 7(Q 2 ) with Q 2 cannot compensate the decrease of a s (Q 2 ) in a 
non-Abelian asymptotically free gauge theory, and the absence of the extra logarithmic dependence in the region 
Q 2 ^ m p would be a serious argument against the asymptotic freedom a s (Q 2 ) <~ l/ln(Q 2 /A 2 ) in favor of the 
assumption [5] that a s (Q 2 ) = const for Q 2 > 2 GeV 2 . 

I am deeply grateful to A.V. Efrcmov for stimulating discussions. 
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NOTES ADDED 

{1} This paper was finished in June 1977, when I was A.V. Efremov's graduate student at Moscow University. The 
original version is reproduced above without changes: just a few typos in formulas were fixed and references are added 
to short notes presented below. They contain updating and clarifying remarks. 

{2} The paper was submitted to Physics Letters B, but was not accepted. After completing my Ph.D. in February 1978, 
we started to write papers containing a detailed description of our approach to pQCD factorization both for inclusive 
and exclusive processes. The first paper [13] presented the outline of the forthcoming cycle [14-17], of which the paper 
on the pion form factor [17] was intended to be the last. It heavily relied on the content of the preceding three papers 
of the cycle. A self-contained paper [18] on factorization for the pion form factor was prepared for the Alushta 1979 
conference. At the conference, I learned from H.D. Politzer about the paper [19] by S. J. Brodsky and G. P. Lepage. 
Their approach was rather different from ours, but they mentioned our papers [13,17]. After their subsequent papers 
[20,21] appeared, we submitted a shortened version [22] of [18] to Physics Letters B. It was published there in 1980. A 
more detailed discussion was published in the same year as the second part of the review in Rivista del Nuovo Cimento 
[23]. Our 1978 paper [17] was eventually also published by Teor. Mat. Fiz. in 1980 (in fact, before the first three papers 
[14-16] of the cycle). 

{3} Nowadays, <^(£,/x 2 ) is usually called the pion distribution amplitude (DA). This terminology was introduced in 
[24] where the pion DA is defined through a k± integral of the light-cone wave function ip(x,k±_). Our definition (10) 
relating the moments of DA to matrix elements of local operators is widely used in modern applications, like higher-order 
calculations, in lattice simulations, etc. 

{4} This simple formula and Eq. (13) were far from being obvious even at the end of 1977, see e.g. Appendix B of [25] 
for an illustration. From my discussions with V. L. Chernyak in the fall of 1977, I learned that he derived an equivalent 
of Eq.(13) well before me, but the relevant publication [26] appeared only in December 1977. 

{5} This form coincides with the asymptotic Q 2 — > oo limit of the pion distribution amplitude (with integration over £ 
reduced to < £ < 1). The explicit statement that the asymptotic form of the pion DA in QCD is given by §/tt(1 — £ 2 ) 
was made in our "advertisement" paper [13]. Implicitly, this result is also contained in the paper by V.L. Chernyak, A.R. 
Zhitnitsky and V.G. Serbo [26] as a statement about the asymptotic limit of their analogue of j(q 2 ). Their derivation 
was outlined in [27]. Our solution of the evolution equation (6) was given in [17]. Under simplifying assumptions (that 
are valid for spin-0 gluons but not in QCD) the solution of the evolution equation was also obtained by Farrar and 
Jackson [28,29]. The (pseudo)scalar gluon model for the pion form factor was studied earlier in our paper [3] (English 
preprint: JINR E2-9717, April 1976). As explained in my review [30], its final result contains the expression amounting 
to the solution of the evolution equation for the pion DA in this model. It is equivalent to the solution given in Refs. 
[28,29]. A very efficient way of solving the DA evolution equations was developed by G. P. Lepage and S. J. Brodsky 
[20]. 
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